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^ Abstract. We study the evolution of a massive black hole pair in a rotationally supported 

I , nuclear disc. The distributions of stars and gas mimic the nuclear region of a gas-rich galaxy 

Q ■ merger remnant. Using high-resolution SPH simulations, we follow the black hole dynam- 

^ ' ics and trace the evolution of the underlying background, until the black holes form a binary. 

C/3 , We find that the gravitational perturbation of the pair creates a core in the disc density pro- 

. file, hence decreasing the gas-dynamical drag. This leads the newly formed binary to stall 

at a separation of ~ 5 pc. In the early phases of the sinking, black holes lose memory of 
their initial orbital eccentricity if they co-rotate with the disc, as rotation of the gaseous 
background promotes circularization of the black hole orbits. Circularization is efficient un- 
til the black holes bind in a binary, though in the latest stages of the simulations a residual 
eccentricity > 0.1 is still present. Black holes are treated as sink particles, allowing for gas 
accretion. We find that accretion strongly depends on the dynamical properties of the black 
holes, and occurs preferentially after circularization. 
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1. Introduction gas-dynamical, gravitational torques excited 

during the merger, when large amounts of 

^ , Collisions of gas rich spiral galaxies may gas is driven into the core of the remnant 

■ trigger starbursts as those observed in lumi- (Kazantzidis et al. 2005; Mayer et al. 2007). 

nous infrared galaxies (LIRGs). A large num- Inside a massive self-gravitating disc, a pu- 
ber of LIRGs hosts a central rotationally sup- tative MBH pair can continue its dynamical 
ported massive (up to 10 "^M©) gaseous disc evolution, and, possibly, can accrete gas, pro- 
extending on scales of ~ 100 pc (Sanders ducing an observable double AGN (Kocsis et 
& Mirabel 1996; Downes & Solomon 1998). al. 2005, Dotti et al. 2006). Here we study the 
These discs may be the end-product of role of a nuclear gaseous disc in driving the or- 
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bital evolution of the MBH binary, and asses 
the possibility of gas accretion into each pair 
member. 

2. Simulation setup 

We follow the dynamics of MBH pairs in 
nuclear discs using numerical simulations 
run with the N-Body/SPH code GADGET 
(Springel, Yoshida & White 2001). 

In our models, two MBHs are placed in 
the plane of a gaseous disc, embedded in a 
larger scale stellar spheroid. The gaseous disc 
is modeled with 2 xlO^ particles, has a total 
mass Moisc - IO^Mq, and follows a Mestel 
surface density profile. The disc is rotation- 
ally supported, and has a radial and vertical 
scale of 100 and 10 pc, respectively. SPH par- 
ticles evolve adiabatically (with a polytrophic 
index = 5/3), and we set the initial internal en- 
ergy density profile so that the Toomre param- 
eter of the disc is > 3 everywhere, preventing 
disc fragmentation and formation of large scale 
over-densities, such as bars and spiral arms. 

The spheroidal component (bulge) is mod- 
eled with 10^ collisionless particles, initially 
distributed as a Plummer sphere, with a core 
radius of 50 pc, and a total mass Msuige - 
6.98Md,sc. 

The two MBHs are equal mass (Mbh - 
4 X 10* Mq). Ml is placed at rest at the cen- 
tre of the circumnuclear disc, while M2 is ini- 
tially orbiting in the plane of the disc on an 
orbit whose eccentricity is ^ 0.7, at a separa- 
tion of 50 pc from Mi . M2 can either be co- or 
counter-rotating with respect to the circumnu- 
clear disc (runs A and B, respectively). 

We allow the gas particles to be accreted 
onto the MBHs if the two following criteria are 
fulfilled: 

• the total energy (kinetic + internal + gravi- 
tational) of the gas particle is lower than 7/10 
of its gravitational energy (all the energies are 
computed with respect to each MBH) 

• the total mass accreted onto a MBH every 
timestep is lower than the M corresponding to 
an Eddington luminosity (LEdd) assuming a ra- 
diative efficiency of 10%. 

The spatial resolution of the hydrodynami- 
cal force in the highest density regions is a; 0.1 



pc. We set the gravitational softening for the 
gaseous particles and the MBHs at the same 
value to prevent numerical errors due to the 
different resolution. With this spatial resolution 
we can resolve the influence radius of Mi 1 
pc), a condition necessary to asses gas accre- 
tion Q 

3. IVIBHs dynamical evolution 

The upper panel of Fig. [U shows the separa- 
tion between the two MBHs as a function of 
time for run A. The two MBHs reach a sep- 
aration of the order of few pc in less than 10 
Myr After an initial fast migration, the shrink- 
ing process becomes inefficient and the binary 
stalls at a separation of a; 5 pc, larger than our 
force resolution. 

The eccentricity evolution of the MBH pair 
is presented in the middle panel of Fig. [T] The 
MBH pair loses memory of its initial eccen- 
tricity in the early phases of the orbital de- 
cay, when the two MBH are at a separation of 
~ 10 pc. Such circularization happens because 
of the different effects that dynamical friction 
exerts on M2 at different orbital phases (see 
Dotti, Colpi & Haardt 2006; Dotti et al. 2007 
for a detailed review). Dynamical friction is ef- 
ficient in reducing the eccentricity down to val- 
ues < 0.1, while, when the binary forms, ec- 
centricity grows again, up to 0.1. 

The stalling of the MBH pair and the resid- 
ual eccentricity are due to the formation of a 
central core in the gaseous disc, as shown in 
Fig. 12] The primary cause of the core forma- 
tion is energy and angular momentum transfer 
from the orbiting MBH to the disc due to the 
dynamical friction. The inner core forms dur- 
ing the early phases of the MBH orbital decay, 
when M2 orbit is still eccentric. We note that 
the interaction between the MBHs and the cir- 
cumnuclear disc in this phase was not fully re- 
solved in Dotti et al. (2007). As a check of our 
initial conditions, we evolved the disc without 
M2 for ~ 10^ yr, and no core formation was ob- 
served. We also ran the same simulation with 
accretion switched off', and we found no differ- 

' The influence radius of depends on the 
phase of its orbit. 
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Fig. 1. Run A. Upper panel: MBH separation as 
a function of time. Middle panel: MBH pair ec- 
centricity evolution. The shaded area refers to the 
time before the first apocentre, when we can not 
evaluate the value of the eccentricity. Lower panel: 
Eddington accretion ratio as a function of time. Red 
and blue lines refer to Ml and M2, respectively. 



10< 



100 



- 








llll|. 
- 




: t=0 










1=2.5 Myr 










1=5 Myr 










; 1=7.5 Myr 
" 1 1 . . 


Mill 






II III" 



I I I I I I ml I I i I ; H 

0.1 1 10 100 



R [pc] 

Fig. 2. Surface density profile of the circumnuclear 
disc at four different times in run A. The solid black, 
blue, red and green solid lines refer to / = 0, 2.5, 5, 
and 7.5 Myr respectively. The vertical dotted lines 
refer to the current position of Mi at the four times. 



ences in the dynamical evolution of the MBH 
pair, as expected due to the small amount of 
gas accreted by the MBHs with respect to the 
total mass of the circumnuclear disc (~ 0.1%). 

For the counter-rotating case (Run B) we 
find that the two MBHs also stall at a separa- 
tion of w 5 pc with a non zero eccentricity as 
illustrated in Fig. |3] The middle panel of Fig. [3] 
shows the evolution of the z-component of the 
orbital angular momentum of M2, normal- 
ized to its initial value. The angular momentum 
(initially negative) grows very efficiently dur- 
ing the first Myr, when M2 is passing through 
the central, high density region of the disc. 
Angular momentum continues to grow mono- 
tonically for the next 3-4 Myrs, then becomes 
positive (Lz w 0.1). That is, M2 starts to move 
on a co-rotating orbit with respect to the disc. 
The dynamical friction process is the ultimate 
cause of this orbital "angular momentum flip". 

4. Accretion processes 

Fig. [T] and |3] allow to compare directly the dy- 
namical properties of the two MBHs to the 



accretion rates in runs A and B respectively. 
The upper panels show the MBHs separation, 
while the lower panels show the Eddington ra- 
tio /Edd = MjM^ii for the central (red lines) 
and orbiting (blue lines) MBHs. In both runs, 
M\ accretes at /edd ~ 0.5, with a slight de- 
crease with time. M2, instead, behaves differ- 
ently. In run A, during the first 7 Myrs, on av- 
erage /Edd ~ 0.4. M2 accretion history can be 
divided in two phases: 

la) for t < 2.5 Myrs the circularization process 
is still efficient. During this phase /sdd ~ 0.3 
on average, showing strong variability; 
2a) for t > 2.5 Myrs, M2 is moving on a quasi- 
circular orbit, and the relative velocity between 
M2 and the gaseous disc is reduced. In this 
phase, /Edd ~ 0.45 on average. 
Larger variability of /Edd can be observed onto 
the counter-rotating M2 in run B. We can still 
distinguish two phases: lb) for t < 3 Myr, M2 
is still counter-rotating (Lz < 0), and /Edd ~ 
0.1, well below /Edd ~ 0.25 obtained averag- 
ing over 5 Myr; 

2b) for f > 3 Myr, M2 accretes at /Edd ~ 0.45 
(on average). 
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Fig. 3. Run B. Upper panel: MBHs separation as a 
function of time. Lower panel: Eddington accretion 
ratio as a function of time. Red and blue lines refer 
to Ml and Mi, respectively. 

5. Conclusions 

Our two high resolution runs show that the 
gas-dynamical interaction between the mas- 
sive circumnuclear disc and MBHs is unable 
to bring the two MBHs at separations of the 
order of the force resolution 0.1 pc) in 5 
Myr. Circularization of the initially eccentric 
orbit of M2 is efficient until the MBHs form 
a binary. In the latest stages of our simula- 
tions 1 eccentricity grows again up to >; 0.1. 
These results are strictly connected to the for- 
mation of a central core in the surface den- 
sity profile of the circumnuclear disc, due to 
the energy and angular momentum exchange 
between the MBH pair and the gaseous par- 
ticles. We stress that we could catch the for- 
mation of the central core (and the stalling of 
the binary) because of the high spatial resolu- 
tion adopted from the beginning of the simu- 
lation, which was increased by a factor x 10 
compared to our previous simulations. It must 
be pointed out that, as discussed in Mayer et 
al. 2007, the fate of the MBH binary could be 
strongly affected by cooling/heating processes, 
not implemented in our runs. 

Thanks to the high spatial resolution of our 
simulations that allows for the code to resolve 



the sphere of influence of the MBHs, we stud- 
ied for the first time the mass accretion rate as 
a function of the dynamical properties of the 
MBHs: we found that variable double nuclear 
activity can be observable for few Myr, when 
the two MBHs orbit with relative separations 
K 10 pc. The accretion rate on counter-rotating 
orbits is more variable, and lower by a factor of 
4-5 when compared to M2 co-rotating with the 
disc. 
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